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Abstract 

We present the results of Monte Carlo simulations of a model of a single maghemite 
ferrimagnetic nanoparticle including radial surface anisotropy distinct from that in 
the core with the aim to clarify what is its role on the magnetization processes at low 
temperatures. The low temperature equilibrium states are analized and compared 
to those of a ferromagnetic particle with the same lattice structure. We have found 
that the formation of hedgehog-like structures due to increased surface anisotropy 
is responsible for a change in the reversal mechanism of the particles. 
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1 Introduction 



The progressive reduction of particle sizes with technological application down 
to the nanometric range results in magnetic behaviours that are predominantly 
dominated by the particle surface spins. The magnetic behaviour of the par- 
ticle surface differs from that corresponding to the core because of the lower 
coordination number and the existence of broken exchange bonds. Among 
other effects, this results in increased anisotropy compared to bulk values 
due to the surface anisotropy with uniaxial normal direction that comes from 
symmetry breaking at the boundaries of the particle [1]. Moreover, depend- 
ing on the character of the exchange interaction between spins in the particle 
and the geometry of their underlying lattice, surface spins may suffer random 
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canting and form a spin-glass state as has been argued in studies of fine par- 
ticles of different ferrimagnetic oxides [2]. These systems display anomalous 
magnetic properties at low temperatures. In particular, they display high clo- 
sure fields with high values of the differential susceptibility and shifts in the 
hysteresis loops after field cooling [3,4,5,6] whose origin is still a matter of 
controversy. While some of these features have been attributed the existence 
of dipolar interactions among the particles [7,8,9,10], there are experimental 
evidences that finite-size and surface effects are crucial in order to understand 
this phenomenology. In this study, we present the results of Monte Carlo (MC) 
simulations of the magnetic properties of individual nanoparticles which aim 
at clarifying what is the specific role played by increased radial anisotropy at 
the surface on the magnetization processes. 



2 Model 



The model considered is an extension to Heisenberg spins of our previous MC 
simulations for an Ising spin model of a 7-Fe 2 3 particle [11,12]. The particles 
considered are spheres of diameter D (in units of the linear unit cell size a). 
The simulations have been run for ferrimagnetic particles with maghemite lat- 
tice structure and compared with those corresponding to ferromagnetic (FM) 
particles with the same underlying lattice. In order to characterize the specific 
contribution of the surface spins to the magnetic properties of the particle, we 
have distinguished two different regions in the particle: the surface formed by 
the outermost unit cells and an internal core of diameter D Cor e — D — 2. 

The magnetic ions are modeled by Heisenberg classical spins with the following 
interaction Hamiltonian: 



N N 

H/k B = ~ E JiA ■ Sj - ]T H ■ S, - ]T [k c {St) 2 + k s (Si ■ n,) 2 } , (1) 

(id) i=l i=l 

which is an extension of our previous study with Ising spins [11]. Here, Si 
are three dimensional classical vector spins placed on a discrete lattice and 
H is the magnetic field, that in the following will be given in temperature 
units as h = fiH/ks, with /i the moment of the magnetic ion. The J^-'s 
are the nearest neigbour exchange constants, which might vary from site to 
site depending on the specific material considered, for the case of maghemite, 
which is a two sublattice (0,T) ferrimagnet, their values are J T t = —21 K, 
Joo — —8.6 K, Jto — —28.1 K [11,4,13]. The last term in Eq. 1 accounts 
for the magnetic anisotropy energy, which has been taken as uniaxial along 
the z-axis for the core spins and along the radial direction for the spins at the 
surface. k c and k s are the core and surface anisotropy constants respectively in 
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temperature units. Although the values of bulk anisotropy constants Kc can 
be obtained indirectly by magnetic measurements, the surface contribution K$ 
to the anisotropy energy is more difficult to evaluate. However, for maghemite 
K c — 4.7 x 10 4 erg/cm 3 [14,15] and K s has been estimated as K s ~ 0.06 
erg/cm 2 from Mossbauer experiments [16,17]. Therefore, in the simulations for 
the maghemite particle, we will vary the core anisotropy values in the range 
kc = 0.01 — IK and those of the surface anisotropy in the range ks = 1 — 100 
K. 

The simulations have been performed using the standard Metropolis algorithm 
[18] for continuous spins. As for the details, let us only mention that, in this 
case, special care has to be taken in the update algorithm with respect to the 
case of Ising spins. In order to sample efficiently the configuration space both 
in the low and high anisotropy limits, we have chosen a sequential combination 
of three kinds of trial angles in a similar way as Hinzke and Nowak [19]. The 
update algorithm consists of a series of three uniform trial steps at random, one 
small trial step inside a cone around the initial spin direction and a reflection 
of the spin direction at the x-y plane to guarantee the correct performance of 
the simulation at high values of the anisotropy. 



3 Results 



In order to gain insight on the equilibrium properties of particles with surface 
anisotropy , we have performed a study of the equilibrium configurations in 
zero magnetic field at low temperatures. For this purpose, the system is ini- 
tially prepared in a high temperature state with spins randomly oriented and 
then it is subsequently cooled in constant temperature steps 5T — 1 down to 
zero temperature. At each temperature we measure the thermal averages of 
the quantities of interest such as the total energy, sublattice magnetizations 
(Mo, M T ), surface M s and core M c contributions to the total magnetization 
MTotai, and also M n = Y^iLi Si ■ hi , the sum of projections of the spins onto 
the local anisotropy axis. 

For a better understanding of the role of surface anisotropy on the zero temper- 
ature configurations, we will first show the results of simulations for a spherical 
particle with D = 3 and kc = 1, and the same structure than maghemite but 
with equal ferromagnetic (FM) interactions (J^ = 1) among the spins instead 
of the real values. In Fig. 1 we show the surface and core contributions to 
thermal dependence of the magnetization during the above-mentioned simu- 
lated annealing procedure for a range of values of ks spanning two orders of 
magnitude. The increase of M z as T is reduced indicates that ordered states 
with FM order start to appear below T ~ 10 K which finally converge towards 
the equilibrium state at zero temperature . As reflected by the M z and M n 
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values attained at T = 0, the increase in ks results in a change in the kind 
of low T order. For ks < 10, quasi-uniform magnetization configurations are 
stable with the core spins pointing along the z axis (M n < 1 and \M Z \ ~ 1 as 
T approaches zero in the right column curves of Fig. 1) and groups of surface 
spins slightly deviated towards the radial direction (M n > 0.5 and \M Z \ ~ 1 
as T approaches zero in the left column curves of Fig. 1). However, as ks is 
further increased from 10, hedgehog-like structures of surface spins start to 
form (note that M n at T = tends to 1 at the surface) inducing the deviation 
of core spins towards the radial direction (indicated by the approach of M n to 
0.5 and of M z to as T decreases for core spins). Snapshots of configurations 
representative of low and high ks values are depicted in Fig. 3 (upper rows) 
in qualitative agreement with simulations performed by Labay et al. for a FM 
particle with sc lattice [20]. 

A similar phenomenon occurs in the case of a ferrimagnetic maghemite spher- 
ical particle with the real antiferromagnetic (AF) exchange constants given in 
the preceding section. However now, due to the dominant AF intersublattice 
coupling, there is a tendency to achieve antiparallel alignment among spins 
in different sublattices that causes a greater degree of disorder at the surface 
with respect to the FM case. Notice also that, even for small ks, M z does not 
approach 1 as T tends to zero as in the FM case, but instead tends to the value 
of magnetization of the noncompensated spins (Mu nc = (No — N^/Nxotai), 
that for a particle of diameter D = 3 is M\j nc = 0.285(0.412) for surface 
(core) spins [11]. Moreover, for the ferrimagnetic particle, there is a change 
in the magnetic order at low T as the k s /k c ratio increases. When surface 
anisotropy is small, the particle orders into a quasi- AF state in which spins in 
each sublattice are almost aligned along the core easy-axis. When increasing 
ks, surface spins start to form throttled structures as those seen for ks in Fig. 
3 that depart the spins from the z-axis. Finally, for ks » k c , hedgehog-like 
configurations are favoured by the dominant radial anisotropy contribution 
(see the configuration for k s = 100). 

In this case, the direct visualization of equilibrium configurations presented in 
Fig. 3 shows that the reduction of the saturation magnetization with particle 
size observed experimentally in different fine particles of ferrimagnetic oxides, 
can be attributed to the random canting of surface spins caused by the com- 
peting antiferromagnetic interaction between sublattices [21,22]. Moreover, as 
the results of our simulations confirm, the degree of disorder at the surface is 
larger than for the FM particle due to the complex interplay between the AF 
intralattice interactions and the local anisotropy easy-axis. 

In order to study the influence of surface anisortropy on the magnetization 
processes, we have also simulated hysteresis loops at zero temperature. The 
loops have been computed by starting from a saturated state achieved after 
application of a high enough field of h — 200 K along the z axis and decreas- 
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ing the field in constant steps 8h = IK, during which the magnetization was 
averaged over 1000 MC steps after thermalization. The results for a spherical 
particle with D = 3 and FM interactions are shown in Fig. 4. For a FM parti- 
cle, the hysteresis loops are dominated by the surface contribution for all the 
values of ks studied as indicated by the non-squaredness of the loops around 
the coercive field. For high values of the surface anisotropy (ks = 50, 100), a 
magnetic field as high as h = 20 K is able to saturate the core but the sur- 
face spins instead point along the radial direction during the magnetization 
process. This is more clearly reflected in the right panels of Fig. 4, where we 
see that for high ks, M n remains close to 1 at the surface during all the rever- 
sal process, while the core spins depart from their easy directions (M n ~ 1) 
dragged by the surface towards the radial direction close to the coercive field, 
where M n ~ 0.5. The hysteresis loops for the FM particle are to be compared 
with those for the ferrimagnetic maghemite particle presented in Fig. 5 for 
ks = 10, 100. First, although the recersal process is still dominated by the 
surface spins, notice that now the loops for the ferrimagnetic particle become 
more elongated and have higher closure fields, ressembling those observed ex- 
perimentally [23]. The increase of ks at the surface results in an increase in the 
coercive field he and a reduction in the high field susceptibility. These features 
are due to the dominance of surface anisotropy over exchange interactions that 
create surface disordered states which become more difficult to reverse by the 
magnetic field. More importantly, there is a change in the magnetization rever- 
sal mechanism when increasing ks from 10 to 100. In the first case (left panels 
in Fig. 5), the core and surface have similar he and closure fields. The particle 
core reverses in a quasi-uniform form with the spins pointing mostly along 
the z-axis (M n ~ 1) except near h c , and with the surface spins following the 
core reversal (with M n < 1 indicating alignment close to the z direction) [24]. 
However, at higher ks (right panels in Fig. 5), surface spins remain close to 
the local radial direction (M n ~ 1) during all the reversal process, driving the 
core spins out of their local easy axis and making their reversal non-uniform 
(0.5 < M n < 1) due to the appearance of the hedgehog-like structures during 
the reversal [24]. 



In conclusion, we have studied the differences in the equilibrium states and 
reversal processes of ferrimagnetic with respect to FM particles, showing in 
both cases that finite-size effects, together with the existence of surface dis- 
order due to increased anisotropy at the surface, are necessary ingredients for 
the low temperature anomalous magnetic properties observed experimentally 
in nanoparticle systems. 
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FIGURE CAPTIONS 



FIG. 1: Spherical (D = 3) particle with FM interactions (Jy = 1): thermal 
dependence of the surface (left column) and core (right column) contributions 
to the projection of the magnetization along the z axis M z and to M n as 
defined in the text during a progressive cooling from a high T at a constant 
rate ST = -1 K. 

FIG. 2: Spherical (D = 3) particle with the real maghemite AF interactions 
quoted in the text: thermal dependence of the surface (left column) and core 
(right column) contributions to the projection of the magnetization along the 
z axis M z and to M n as defined in the text during a progressive cooling from 
a high T at a constant rate ST = — 1 K. 

FIG. 3: Snapshots of the equilibrium configurations attained after cooling from 
a high temperature disordered state for FM (upper row) and ferrimagnetic 
maghemite (lower row) spherical particles with D = 3 and maghemite lattice 
structure. 

FIG. 4: Hysteresis loops for a spherical particle of diameter D = 3 with FM 
interactions. The core anisotropy constant is kc — 1, the results for several 
values of the surface anistropy constant are displayed k s = 1, 5, 10, 20, 50, 100. 
The magnetization component along the applied field direction is shown in 
lower right panel, the core contribution is displayed in the upper right panel. 
Right panels show the suface and core contributions to M n as defined in the 
text. 

FIG. 5: Hysteresis loops for a spherical particle of diameter D = 3 with the 
real values of exchange constants of maghemite. The main panels show the 
total magnetization along the field direction and the corresponding surface 
(dashed lines) and core (circles) contributions. The upper subpanels show the 
surface (circles) and core (continuous lines) contributions to M n . Left panel: 
kg = 10, right panel: ks = 100. 
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Fig. 3. 
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Fig. 4. 



12 



This figure "Figure3.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/0307584vl 



